function c2 = eq_c2_1(T,W)
% c2 = eq_c2_1(T,W)
% Fitted c2(T,W) as in Dong's thesis
%
% Reference:
% [1] Dong's thesis

k_c = 1.0274e-13;
a_c = 0.1669;
E_c = 1.9e-3; % Unit = EV
k_B = 8.617343e-5; % Boltzmann's constant, Unit = ev/K

c2 = k_c*exp(k_B*T/E_c)*W^(-a_c)/(100*100);

